Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 



(MN WF$t style file v2.2) 



CMB POLARIZATION DUE TO SCATTERING IN 
CLUSTERS 



o 
o 

(N 
(N 



> 
oo 

in 

(N 

o 

O 

Or 

6 



M. Shimon 1 , 2 , Y. Rephaeli 1 , 2 , B. W. O'Shea 2 , 3 , 4 and M. L. Norman 2 

1. School of Physics and Astronomy, Tel Aviv University, Tel Aviv, 69978, Israel 

2. Center for Astrophysics and Space Sciences, University of California, San Diego, 9500 Gilraan Drive, La Jolla, CA, 92093-0424 

3. Department of Physics, University of Illinois at Urbana- Champaign, 1110 West Green Street, Urbana, IL 61801-3080 

4- Theoretical Astrophysics Group, Mail Stop B227, P.O. Box 1663, Los Alamos National Laboratory, Los Alamos, NM 87545 



ABSTRACT 

Scattering of the cosmic microwave background (CMB) in clusters of galaxies 
polarizes the radiation. We explore several polarization components which have their 
origin in the kinematic quadrupole moments induced by the motion of the scattering 
electrons, either directed or random. Polarization levels and patterns are determined 
in a cluster simulated by the hydrodynamical Enzo code. We find that polarization 
signals can be as high as ~ 1 /xK, a level that may be detectable by upcoming CMB 
experiments. 
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1 INTRODUCTION 

Scattering of the CMB in clusters of galaxies imprints 
spectral changes on the Planck spectrum - the Sunyaev- 
Zeldovich (S-Z) effect - that constitute major cluster and 
cosmological probes. Cluster properties, such as intracluster 
(IC) gas density, temperature, gas mass, total mass, and pe- 
culiar velocities, as well as values of the basic cosmological 
parameters, can be determined from measurements of the 
effect in a (sufficiently) large sample of clusters (as reviewed 
by Rephaeli 1995, Birkinshaw 1999, and Carlstrom et al. 
2002). While this wealth of astrophysical and cosmological 
information can be extracted largely from measurements of 
the total intensity change due to Compton scattering, the 
polarization state is also of considerable interest. Various 
polarization components have been identified in the original 
work of Sunyaev & Zeldovich (1980), and further elaborated 
upon by Sazonov & Sunyaev (1999). Measurements of polar- 
ization patterns in clusters may yield additional information 
on the dynamical state of the cluster, the cluster morphol- 
ogy, and its velocity component transverse to the line of 
sight (los). 

Linear polarization is produced by Compton scatter- 
ing when the incident radiation either has a quadrupole 
moment, or acquires it during the (first stage of the) scat- 
tering process. The level of S-Z polarization is low since it 
is proportional to the product of the quadrupole moment 
and the Compton optical depth, r. Cluster polarization sig- 
nals are predicted to be comparable to or below ~ 1 fiK. 
The detection of polarized CMB signals at this level seems 



now feasible, and indeed anticipated by several experiments 
(e.g. Bowden et al. 2004). The projected resolution and 
sensitivities motivate investigation of the polarization levels 
and patterns that are expected in clusters; this is the ba- 
sic objective of the work described here. Weak polarization 
components that are of interest, but ones that are not ex- 
plored here, are produced when the radiation is scattered 
during aspherical collapse of protoclusters, scattering in a 
moving cluster in which the radiation develops anisotropy 
due to lensing (Gibilisco 1997), or polarization produced by 
IC magnetic fields (Ohno et al. 2003) and relativistic mag- 
netized plasma (Cooray, Melchiorri & Silk 2002). 

The dynamical state of clusters undergoing mergers of 
subclusters cannot be adequately represented by simple an- 
alytical models. Since the polarization levels and patterns 
are quite sensitive to the non-uniform distributions of the 
gas and total mass, peculiar and internal velocities, and the 
gas temperature profile, a realistic characterization of the 
polarization properties can only be based on the results of 
hydrodynamic simulations. In the first stage of a simulation- 
based study, we present here results for the S-Z polarization 
properties in clusters simulated by the hydrodynamic Enzo 
code. We illustrate these for a simulated 1.4 x 10 15 M@ clus- 
ter. 

In the next section we summarize the cluster simula- 
tion work whose results are employed here. Section 3 briefly 
describes the various polarization components induced by 
Compton scattering in clusters. In section 4 we describe our 
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simulation-based results, and end with a discussion in sec- 
tion 5. 



2 SIMULATED CLUSTERS 

In this paper we use the data obtained from Enzo, an Eule- 
rian adaptive mesh cosmology code (O'Shea et al. 2004 and 
references therein). This code solves dark matter N-body 
dynamics with the particle mesh technique. The Poisson 
equation is solved using a combination of fast Fourier trans- 
forms on the periodic root grid and multigrid techniques 
on non-periodic subgrids. Euler's hydrodynamic equations 
are solved using a version of the piecewise parabolic method 
(PPM) algorithm of Woodward & Colella (1984) which has 
been modified for high-Mach flows such as those found 
in galaxy cluster simulations (Bryan et al. 1995). Block- 
structured adaptive mesh refinement using the technique 
of Berger & Colella (1989) is used to follow the evolution 
of objects of interest. Enzo allows arbitrary integer ratios of 
parent and child grid resolutions and mesh refinement based 
on a variety of criteria, including baryon and dark matter 
overdensity or slope, the existence of shocks, Jeans length 
resolution, and cell cooling time. The code is parallelized 
using the MPI message-passing library . Additional physics 
packages are available (but not used for the simulations in 
this paper) and are described in O'Shea et al. (2004). The 
Enzo code is publicly available on the world wide web 2 . 

The simulation was initialized as follows: In order to 
find a galaxy cluster of appropriate size we initialized a 
simulation using only dark matter in a volume of space 
256 h" 1 Mpc on a side at z = 30 (when density fluctua- 
tions are still linear on scales relevant here). The simula- 
tion used a 256 3 root grid with 256 3 dark matter particles 
(M DM = 1.19 x 10 11 M Q ). The standard ACDM cosmolog- 
ical model was used with Q to t = 1, Ha = 0.7, 57m = 0.3, 
Hq = 70 km/sec/Mpc, as = 0.9, and n = 1. Adaptive mesh 
refinement was turned on with an additional four levels of 
refinement (doubling the spatial resolution at each level) and 
the simulation was run to z = 0. At this point we stopped the 
simulation and used the Hop halo-finding algorithm (Eisen- 
stein & Hut 1998) to find the most massive dark matter halo 
in the simulation. This halo has a mass of 1.4 x 10 15 M© at 

2 = 0. 

The simulation was then re-initialized at z — 30 with 
both dark matter and baryons (Qt — 0.04, SIcdm = 0.26) 
with a 128 3 top grid and 2 static nested grids covering the 
Lagrangian volume where the cluster forms, giving an effec- 
tive root grid resolution of 512 3 cells (0.5 Mpc/h comoving) 
and dark matter particles of mass 1.49 x 10 10 M Q . Adaptive 
mesh refinement is allowed in the region where the galaxy 
cluster forms with a total of 7 levels of refinement beyond 
the root grid, for a maximum spatial resolution of 15.625 
h^ 1 kpc ( comoving). This simulation was then evolved to 
z = 0.06, following the evolution of the dark matter and 
using adiabatic gas dynamics. The galaxy cluster of inter- 
est was then relocated with the Hop algorithm at z — 0.06, 
0.25, 0.5 and 1.0, and data cubes with 256 3 cells containing 
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Pdm, Pgas, T e and the three baryon velocity components 
were extracted at these redshifts. The total cube length is 
4 Mpc/h, corresponding to a spatial resolution of 15.625 
kpc/h (comoving). Unless otherwise noted, the results in 
the following sections are for the full 256 3 data cubes, and 
are based on the simulated cluster at z = 0.5. The baryon 
density and velocity field are shown in Figure 1 and the clus- 
ter projected temperature and optical depth are depicted in 
Figure 2. 



3 CLUSTER POLARIZATION COMPONENTS 

The main process that polarizes the CMB in clusters is 
Compton scattering of the radiation by IC electrons. Polar- 
ization levels depend on the cluster optical depth to Comp- 
ton scattering, and either on the cluster velocity or gas tem- 
perature. These are typically in the 10 — 100 nK range, but 
may reach values as high as ~ ^iK. We describe the main 
polarization components, and illustrate their predicted pat- 
terns in our fiducial cluster. As we have noted, the CMB can 
be polarized by Compton scattering if the radiation multi- 
pole expansion has a finite quadrupole component in the 
electron rest frame. This can be brought about in several 
ways, as described in the following subsections. The level of 
linear polarization and its orientation are determined by the 
two Stokes parameters 
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where dl is a length element along the photon path, n e is the 
electron number density, <tt is the Thomson cross section, 
and a and ip define the relative directions of the incoming 
and outgoing photons. The average electric field describes 
the polarization plane with a direction given by 
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When the incident radiation is expanded in spherical har- 



I(a,ip) = ^2 Ii m Yi m (a,ip), 



(3) 



it becomes evident (from the orthogonality conditions of 
the spherical harmonics) that only the quadrupole moment 
terms proportional to 72,2 and 72,-2 contribute to Q and U 
in Equation (1). 

In the following subsections we describe the polariza- 
tion produced by Compton scattering when a quadrupole 
moment is induced by electron motions. These 'kinematic' 
polarization components arise due to either bulk (peculiar) 
or thermal (random) motion of the electrons. Since we are 
mainly interested in a more realistic assessment of the vari- 
ous components based on the results of hydrodynamic sim- 
ulation (rather than idealized analytic models), previously 
derived theoretical expressions are only briefly summarized 
here. 
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3.1 POLARIZATION INDUCED BY THE 
GLOBAL CMB QUADRUPOLE 

The CMB global quadrupole moment originates in the 
Sachs- Wolfe (SW) effect at the surface of last scattering, and 
later in the evolving gravitational fields of density inhomo- 
geneities - the integrated Sachs- Wolfe (ISW) effect. Scatter- 
ing by IC electrons polarizes the radiation at a level that is 
(linearly) proportional to the rms of the primary quadrupole 
moment and the cluster optical depth to Compton scatter- 
ing, r (Sazonov & Sunyaev 1999). Using the WMAP nor- 
malization of the CMB quadrupole moment (Bennett et 
al. 2003) the maximal polarized signal is expected to be 
~ 2.6r nK, and its all-sky average is ~ 60% of this value 
(Sazonov & Sunyaev 1999). The polarization pattern has 
broad peaks extending over an appreciable fraction of the 
sky. The minima are narrow and correspond to the change 
of sign of the Stokes parameters. It has been noted that 
the dependence of this polarization component on the CMB 
quadrupole moment could possibly be exploited as a means 
to reduce cosmic variance (Kamionkowski & Loeb 1997), and 
as a probe of dark energy models through the redshift evolu- 
tion of the quadrupole (see, however, Portsmouth 2004). We 
note that the polarization due to the CMB quadrupole can 
be distinguished from the other cluster polarization compo- 
nents by virtue of its large-scale distribution reflecting that 
of the primary CMB quadrupole (e.g. Baumann & Cooray 
2003), and the fact that it is independent of frequency (when 
expressed in temperature units). This polarization compo- 
nent, which was recently studied in the context of N-body 
simulations by Amblard & White (2004), will not be dis- 
cussed further here. 



asymmetry introduced by the direction of the cluster mo- 
tion. 

To calculate the polarization components, we first note 
that the expressions for the Q and U Stokes parameters are 
specified in Equation (1) in terms of the relative directions 
of the ingoing and outgoing photons. Due to the obvious 
dependence of the angles and frequencies on the relative di- 
rections of the electron and photon motions, the expressions 
for Q and U are transformed to a frame whose Z axis coin- 
cides with the direction of the electron velocity, the axis with 
respect to which the incoming and outgoing photon direc- 
tions are defined (in accord with the choice made by Chan- 
drasekhar 1950). Since the Doppler effect depends only on 
the angle between the photon wave vector and the electron 
velocity, the calculation of the temperature anisotropy and 
polarization is easier if we assume the incident radiation is 
isotropic in the CMB frame (which is indeed the case to the 
level of accuracy relevant to this work). With the polarized 
cross section for Compton scattering (in the Thomson limit) 
written in terms of angles measured with respect to the elec- 
tron velocity (after averaging over the azimuthal angles), 

£ = ^(i-/4!)fl.04), W 

where P2(Mo) i s the second Legendre polynomial, the expres- 
sion for the Q parameter of the scattered radiation is 

3 2 f 1 

Q{v) = gf(l - A*o) J P2(n'o)I(n'o)dn' . (5) 

Here, f/ = cos0' o , /no = cos#o, are the angle cosines be- 
tween the electron velocity and the incoming and outgoing 
photons, respectively. 



3.2 POLARIZATION DUE TO INTRINSIC 
QUADRUPOLE MOMENTS 

Scattering of the CMB in a cluster results in local anisotropy 
due to the different pathlengths of photons arriving from 
various directions to a given point. The optical depth effec- 
tively develops a dipole moment, and AT/T would appear 
anisotropic, except in the center of a spherically symmet- 
ric cluster, where electrons effectively 'see' almost isotropic 
(singly) scattered radiation. The local anisotropy of AT/T 
resulting from single scatterings is due to either thermal or 
bulk electron motions, constituting the well known thermal 
and kinematic components of the S-Z effect. This anisotropy 
provides the requisite quadrupole moment; second scatter- 
ings then polarize the radiation. 

Two polarization components are directly induced by 
scattering in a moving cluster. The first is linear in the clus- 
ter velocity component transverse to the los, v t = c(3 t , but 
quadratic in r; the second is linear in r but quadratic in fit- 
A third polarization component arises from second scatter- 
ings (and, therefore, is oc r ) after local anisotropy is pro- 
duced by the thermal S-Z effect. The spatial patterns of the 
various polarization components can be readily determined 
only when the morphology of the cluster and its IC gas are 
isotropic. The polarization patterns arising from scattering 
off thermal electrons are isotropic (Sazonov & Sunyaev 1999) 
in a spherical cluster, while the corresponding patterns of 
the kinematic components are clearly anisotropic due to the 



3.2.1 The t fit and,Tf3 t Components 

The degree of polarization induced by the kinematic S-Z 
component was determined by Sunyaev & Zeldovich (1980) 
in the simple case of uniform gas density, 
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where x = hv/kT is the photon dimensionless frequency, 
with T denoting the incident CMB temperature. A more 
complete calculation of this and the other polarization com- 
ponents was formulated by Sazonov & Sunyaev (1999). 
Viewed along a direction h = (9,<fi), the temperature 
anisotropy at a point (X, Y, Z), AT(X, Y, Z, 8, 4>), generates 
polarization upon second scatterings. The Stokes parame- 
ters are calculated from Equation (1), 
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where now the Z direction is chosen along the los, and To is 
the incident intensity. The intensity change resulting from 
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first scatterings is 



AI{X,Y,Z,0,</>)/I o = a T / dl(X',Y',Z',8,<t>) 



xn e (X', Y', Z', 6, 4>)h ■ p{X', Y , Z') 
xe x 



1 



(8) 



and the optical depth thru the point (X, Y, Z) in the direc- 
tion (6, <f>) is 



(9) 



T(X,Y,Z,e,ct>) = o T J n e (X',Y',Z')dl(X',Y',Z',e, 

Q(X,Y) and U(X,Y) fully describe the 2-D polarization 
field. 

The second kinematic polarization component is oc r/3 2 ; 
this component is generated by virtue of the fact that the ra- 
diation appears anisotropic in the electron frame. Expanding 
the apparent angular distribution of the radiation in Legen- 
dre polynomials, and keeping terms up to f3 2 , the quadrupole 
moment is determined. The polarization of the singly scat- 
tered radiation is then calculated (Sazonov & Sunyaev 1999) 
by using Equation (5). (Note that the transformation back 
to the CMB frame does not change the polarization since 
this introduces additional terms which are higher order in 
(5.) The level of this polarization component (which was orig- 
inally determined by Sunyaev & Zeldovich 1980) is 

In our chosen frame of reference the Stokes parameter U van- 
ishes due to azimuthal symmetry, so the total polarization 
amplitude, P, is equal to Q and the polarization is orthog- 
onal to fit- Relativistic corrections (to the non-relativistic 
expression of Sunyaev & Zeldovich 1980) were calculated 
by Challinor, Ford & Lasenby (2000) and Itoh et al. (2000). 
These corrections generally amount to a ~ 10% reduction in 
Q. In our calculations we use their analytic approximation 
(to the exact relativistic calculation) 

Q = 10(c* - l) 2 TPt 

^F + e(3F-2 (2F 2 + G 2 ) + ^F (f' 2 + 2G 2 



+ 



--F+ - (2F 2 +G 2 



(11) 



where F = x coth | , G = x/ sinh | , Q is the gas dimension- 
less temperature, and c(3 r is the radial component of the 
cluster velocity. 



3.3 POLARIZATION DUE TO SCATTERING 
BY THERMAL ELECTRONS 

Analogously to the r 2 /3 component discussed above, dou- 
ble scattering off electrons moving with random thermal ve- 
locities can induce polarization whose level is proportional 
to r 2 0. The anisotropy introduced by single scatterings is 
the usual thermal component of the S-Z effect with inten- 
sity change A7 t . We use an analytic approximation to the 
exact relativistic calculation for AI t that was obtained by 
Itoh, Kohyama & Nozawa (1998), and Shimon & Rephaeli 



(2004). Using the expression for AI t (from the latter pa- 
per) in Equation (7), the polarization can be calculated as 
discussed above 



A/ t (x,y,z,M) 



utIox 



n e (r)9(r)'dZ. (12) 



In Equation (12) the integration is along the photon path 
prior to the second scattering, fi(x) are spectral functions 
of x that are specified by Shimon & Rephaeli (2004). 



4 RESULTS 

We have calculated the polarization patterns of the kine- 
matic and thermal components using a simulation of a mas- 
sive cluster at z = 0.50, with M c i — 1.4 x 10 15 Mq and veloc- 
ity dispersion of » 1000 km/s. The baryon density and the 
velocity field in the relevant region of the simulated cluster 
are shown in Figure 1. Polarization levels in all the figures 
are specified in terms of the equivalent (absolute) tempera- 
ture. The figures show results in the sky plane with the los 
along the Z-axis. 

The calculations of the double scattering components 
(proportional to r 2 and t 2 (3) are more challenging than 
that of the single scattering case. The calculations generally 
require ray tracing of doubly scattered photons. In Cartesian 
coordinates, Equations (7) & (12) can be written as 

Q(X,Y) = M J J dZd 3 r'n e (rK(r') 



(X- X') 2 -(Y -Y' 



U(X,Y) 
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where the first integration (over Z) is along the los. The 
electron density n e was calculated from the simulated p g 
taking a molecular weight of fi — 0.63. 

Calculation of these 4-D integrals over the simulation 
256 3 data cells is computationally intensive, motivating an 
attempt to reduce the level of computation. For example, 
an approximation (employed by Lavaux et al. 2004) is to 
describe the first scattering by representing the gas in terms 
of a best-fit /3 - profile. This significantly reduces the extent 
of the calculation. From Figures 1 & 2, it is evident that our 
cluster cannot be fitted by a (3 profile, isothermal or not; 
fitting merging clusters with subclumps and rich substruc- 
ture with a p profile is clearly inappropriate. Fortunately, 
Equation (13) can be cast in the form 



P+AX, Y) 



= a 2 T J 
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(14) 



where P + and P x are the Stokes Q & U parameters, respec- 
tively, and 



= /+,x *9 



(15) 
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is a convolution of the functions / +jX with g, where 
X 2 -Y 2 _ 2XY 

g = ^n e (r)/,( 2: )e(r) 1 . (16) 

i 

This representation enables us to employ 3-D FFT calcu- 
lations that are considerably faster than the direct integra- 
tions of Equations (13). However, the full computation is 
still quite extensive; we therefore focus here on the core of 
the simulation volume. 

Our results are summarized in Figure 3. Unlike the case 
of spherical (Sazonov & Sunyaev 1999) or quasi-spherical 
(Lavaux et al. 2004) clusters, we obtain a more intricate 
polarization pattern and the quadrupole structure, which is 
absent in the spherically symmetric case, is evident. This 
pattern reflects the non-uniform gas distribution and an ap- 
preciable degree of non-sphericity of the cluster, which sub- 
stantially affect primarily the double scattering polarization 
signals, essentially due to the fact that the second scatter- 
ing is more likely to occur far from the cluster center (see 
also Sazonov & Sunyaev 1999, Lavaux et al. 2004). The 
level of polarization at the 545 GHz Planck frequency band, 
at which the beam is 4.5' FWHM, can reach ~ 80 nK, but 
when convolved with a narrow beam 1' FWHM profile, its 
peak value can be as high as ~ 0.7fiK. This level definitely 
grazes the threshold of next generation experiments. Note, 
however, that we have considered a rich cluster with max- 
imal optical depth of To £s 0.01 (Figure 2), but since this 
component is proportional to ~ r 2 , typical values could be 
a factor ~ 4 lower or higher than the specific values quoted 
here. Additional uncertainty stems from the scaling of the 
cluster temperature with mass, and from the evolution of 
the temperature with redshift. Finally, this component was 
calculated at 545 GHz; measurements at lower frequencies 
will result in lower levels of polarization (Equation 13). 

Equations (7)- (9) for the double scattering kinematic 
component (oc r 2 /3) can be written in the same form as 
Equation (14) , but with 

3 

h+, x = /+,x *g = 5^/i+, x *9i (17) 

i=l 

where 

X^-Y 2 ) 2X l XY 
fi+ = -g , fix = — -g — , Si = n(r)/3i(r) (18) 

with Xi (i — 1, 2, 3) the components of the vector (X, Y, Z). 
The results for the r 2 /3 component (for the central region 
of the simulation) are summarized in Figure 4. The level of 
polarization in absolute temperature units is ~ 1 — 2 orders 
of magnitudes smaller than the oc r 2 component. Con- 
volving with the Planck HFI and a gaussian beam profile 
with FWHM of 1', the peak values are 3 nK and 25nK, 
respectively, comparable to the level of the cosmological 
quadrupole polarization (section 3.1). As previously noted, 
both these components are frequency-independent, a fact 
that has obvious implications on the feasibility of separating 
these components out when the measurements will eventu- 
ally reach the requisite sensitivity. 

The polarization at 545 GHz due to the Doppler 
anisotropy, oc r/3 2 , is plotted in Figure 5. At this high fre- 
quency the polarization can be as high as ~ 18 nK when 



convolved with the Planck HFI. Convolution with a narrow 
gaussian with 1' FWHM beam yields higher (~ 0.1/iif) lev- 
els. Owing to its spectral dependence, this component can 
- in principle - be distinguished from both the oc t 2 © and 
oc t 2 (] polarization by multi-frequency measurements. More- 
over, due to the different dependence on the optical depth, 
the two kinematic effects (oc t 2 /3 and oc r/3 2 ) can possibly be 
separated out through their different morphologies (compare 
Figures 4 & 5). 



5 DISCUSSION 

CMB observational capabilities have dramatically improved 
since the first measurement of the primary anisotropy by 
COBE/DMR. With projected sensitivities around 1 fjK, and 
sub-arcminute angular resolution, measurement of the dom- 
inant polarization signals in clusters by upcoming experi- 
ments seems feasible in the (not too distant) future. 

The coupling of the CMB to IC gas by Compton scat- 
tering leads to various polarization components whose mag- 
nitudes and spatial patterns could provide significant infor- 
mation on the cluster motion and the dynamical state of IC 
gas. The feasibility of identifying a given polarization com- 
ponent (obviously) depends primarily on its magnitude, but 
this will be aided also by its particular spectral and spatial 
characteristics. The fact that some of the spatial patterns 
(such as that of the thermal S-Z polarization) produce no 
net signal (or substantially reduced signal in aspherical clus- 
ter) when integrated over a sufficiently large region of the 
cluster surface, can be exploited in optimizing the observa- 
tional strategy for measurement of a specific polarization 
component. 

In this work we have used cluster simulations to assess 
the amplitude of three CMB polarization components. As 
expected, most polarization signals induced by the interac- 
tion of CMB photons and IC gas are substantially smaller 
than 1 /iK. An exception is the polarization induced by the 
double scattering thermal effect. For this component, con- 
volving our simulation data with a FWHM 1' beam, we find 
that the effective polarization signal is about ~ 1/j.K. The 
polarization induced by double scattering off the thermal gas 
was the focus of the work by Lavaux et al. (2004). We have 
implemented a direct calculation scheme that yields results 
at the cluster core which qualitatively differ from those of the 
latter authors. By assuming isothermal /3 profile to describe 
first scatterings, Lavaux et al. (2004) partially underesti- 
mated the effect of dynamics on producing the quadrupole 
in the first scattering; this may be a reason for the qual- 
itative discrepancy between our respective results. By ex- 
ploiting the properties of convolution integrals we are able 
to considerably facilitate the calculation of the components 
oc t 2 0, oc t 2 /3, essentially avoiding the need for simplifying 
approximations. The oc r 2 /3 and oc r/3 2 components are an 
order of magnitude weaker, too small to be detected towards 
individual clusters in near future experiments. 

This is the first stage of a comprehensive study of CMB 
polarization in clusters based on a realistic description of 
cluster structure that is afforded by hydrodynamic simula- 
tions. The results presented here already point to apprecia- 
ble discrepancy with what is obtained when IC gas is mod- 
eled by a simple spherical distribution. We plan to extend 
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the work presented here to a sample of simulated clusters 
with a range of masses and redshifts. The polarization sig- 
nals in each of the sampled clusters together, supplemented 
by the expected confusing polarized signals, will then be an- 
alyzed in order to determine optimal strategies for the mea- 
surement of cluster-induced polarization components and 
the extraction of cluster properties. 
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Figure 1. The velocity field is shown on the background of the 
baryon mass density; arrow length is linearly proportional to ve- 
locity magnitude. Color density scale is in units of g/cm 3 . 
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Figure 2. The gas temperature (top figure, with color scale in 
K) and optical depth (bottom) 



CMB POLARIZATION DUE TO SCATTERING IN CLUSTERS 7 



■nponsntat 545 GHz i Planck H Fl J 

wwCa vw v \ \ ^ v W 
WWW w 
\\ \ w\\w\ 
w\ww 
<fc\ vww 

S£\ \ \\ w 

\\\\s w 

\ \ w \ s 




t B corTpon?nl a1 545 GH: iP/VHM 1 sicminbgaml 




t'P polarization ccnponent [Planck H Fl J 




S/"*"SSS/-S /S S / / , 

f-v J- S f s s s s s s s > y w f / s J f / -4 * 
//./ ^ /////■ /-V i 



t"P component i"FvVHIj! 1 aicmln baamj 




Figure 3. The r 2 polarization component convolved with the 
Planck HFI beam (top) and with a FWHM 1' beam (bottom). 
Color scale is in /zK. 



Figure 4. The r 2 /3 polarization component convolved with the 
Planck HFI beam (top) and with a FWHM 1' beam (bottom). 
Color scale is in fiK. 
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Figure 5. The r/3 2 polarization component at 545 GHz convolved 
with the Planck HFI beam (top) and with a FWHM 1' beam 
(bottom). Color scale is in /iK. 



